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A well known description of superradiance from pointlike collections of many atoms involves 
the dissipative motion of a large spin. The pertinent "superradiance master equation" allows for a 
formally exact solution which we subject to a semiclassical evaluation. The clue is a saddle-point 
approximation for an inverse Laplace transform. All previous approximate treatments, disparate as 
they may appear, are encompassed in our systematic formulation. A byproduct is a hitherto un- 
known rigorous relation between coherences and probabilities. Our results allow for generalizations 
to spin dynamics with chaos in the classical limit. 
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I. INTRODUCTION 



Dissipative motion of large spins was first seen in experiments on superradiance or superfluorescence (For extensive 
reviews see Refs. @^]), after being proposed a lot earlier by Dicke ||. The so called superradiance master equation 
proposed in Q|| has since become a standard tool for describing the collective dynamics of identical superradiating 
atoms in the small-sample limit. Formally speaking, it provides a quantum treatment of a large spin with conserved 
square, J 2 = j(j + 1), with the quantum number j capable of taking on positive half integer or integer values up to 
half the number of atoms N. The origin of such an angular momentum lies in the familiar formal equivalence of a 
single two-level atom to a spin-i. In (semi) classical parlance, the spin in question is called the Bloch vector whose 
z-component measures the energy stored in atomic excitation while the transverse components are related to the 
dipole element responsible for the atomic transition. More or less everything worth knowing about the superradiance 
master equation in relation to the numerous superfluorescence experiments has been worked out more than a decade 
ago. 

When we pick up the thread now our motivation is not to better explain anything previously observed, but rather 
the expectation of new experiments involving dissipative motion of large spins constituted by many identical two-level 
atoms, albeit motions that would have a chaotic classical limit and display quantum manifestations of chaos when 
the spin quantum number j is of the order of several hundreds or thousands. When beginning to look into such 
dynamics Q we found, somewhat to our surprise, that previous treatments of the superradiance master equation were 
so directly geared to the specifics of superradiant pulses as transient events that new questions do indeed require some 
new theoretical work. In particular, the semiclassical limit of large j deserves systematic attention and turns out to 
harbor one or the other surprise which we begin to uncover in the present paper. 

The large- j limit can be approached through the rigorous solution of the master equation which was known from the 
very beginning Q, and we shall actually follow that path here. Strangely enough, up to now that rigorous solution 
has mostly been looked upon as a curiosity rather than a useful starting point of analytic work; even numerical 
evaluations were disfavored against routines for solving coupled differential equations for density matrix elements in 
some representation. 

We propose to show that the large- j limit is very conveniently accessed by subjecting the rigorous Laplace trans- 
formed density matrix to a saddle-point evaluation of the inverse Laplace transformation. More specifically, we carry 
out this program in the eigenrepresentation of J z and J 2 for the density matrix (jm\ p(t)\jm') and the propagator 
relating that density matrix to its initial form (jm|p(0)|jm'). The saddle-point result turns out reliable provided that 
not only j is large but also the difference between the initial and final eigenvalues of J z , i.e. \m — m'\ 1. That 
restriction unfortunately affects the propagator at early times while most of the probability still resides in levels m 
close to the initial m' . We therefore establish an independent early-time propagator, show its agreement with the 
saddle-point version in a certain time span and finally combine the two to an explicit expression of uniform validity. 

Our uniform propagator turns out to systematically encompass previous asymptotic results. Among these is, 
trivially, the fully classical behavior arising in the limit j ' ^ oo as long as the initial state is not too close to the state 
of full initial excitation m = j which in the classical limit is an infinitly long-lived state of marginal equilibrium. The 
classical behavior in question is that of an overdamped pendulum. The pertinent equation of motion for the so-called 
Bloch angle 9 (defined through cos 9 = lim :) ^ 00 ( J z (r))/j) reads, with r denoting a suitably scaled time, -^9 = sin9; 
the the well known solution is 

tan^=e^tan^. (1.1) 
2 2 V ' 

Furthermore, we recover the random-jitter picture first suggested in and the ensuing distribution of delay times 
as well as the scaling results for time dependent expectation values of products of the observables J x , J y , J z obtained 
by somewhat hit-and-run methods in j9|,|To[| . 

An interesting byproduct of our investigation is an exact relation between diagonal and offdiagonal elements of 
the density matrix in the jm-basis, which to the best of our knowledge has previously gone unnoticed. One may 
thus confine all work towards solving the master equation to the probabilities (jm\p(t)\jm) and eventually obtain the 
coherences (jm\p(t)\jm') through the relation in question. 

A subsequent paper will deal with the large- j limit with the WKB method. 

II. MASTER EQUATION AND DISSIPATIVE PROPAGATOR 

The two states of an atom resonantly coupled to a mode of the electromagnetic field may be thought of as the 
states of a spin-i, and all observables of the effective two- level atom can be represented as linear combinations of 
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unity and the three spin operators J x , J y , J z . In particular, the energy may be associated with J z and the other two 
spin operators with the atomic dipole moment. If N such atoms, all identical, couple collectively to the electric field 
E one has an interaction Hamiltonian cx — J X E where J x = Ylu=i J x is the sum of all single-atom contributions; 
similarly, one has a global atomic energy cx 

Jz = E,Li J z- The collective spin operators obey the familiar angular- 
momentum commutation relations [J x , J y ] = i,J z etc. The Hilbert space for the N atoms is 2 N dimensional but falls 
into subspaces not connected by the collective observables J^; each subspace has fixed J 2 = j(j + 1) with nonnegative 
integer or half-integer j not exceeding N/2. The (2j + 1) states in the jth subspace are conveniently taken as the 
eigenstates \jm) of J z with eigenvalues m — —j, —j + 1, —j + 2, . . . , j. The highest energy may be associated with 
m = j whereupon the ground state has m = —j. In particular, the subspace with J = consists of N + 1 states 
which are all totally symmetric in all atoms; that space may be singled out experimentally by preparing all atoms in 
their lower state. 

In the superradiance experiments of Ref. [[llj a single mode of the electromagnetic field within a resonator was cou- 
pled to N two-level atoms such that the dynamics was that of the so-called Jaynes-Cummings model, with dissipation 
included to account for field losses from the resonator. In the limit of overdamped Rabi oscillations the field mode 
can be eliminated adiabatically. A master equation for the atomic density operator thus results of which we shall 
consider the low-temperature version, thus forbidding the atoms to pick up thermal photons from the environment, 

j t P = K{[J-,pJ+] + [J-P,J+]}\ (2.1) 

here J± = J x ± iJ v are the familiar raising and lowering operators and k measures the rate of photon loss from the 
cavity. 

In the basis set \jm) we obtain from (p.Ij ) a set of equations for the elements p mi m 2 — (j m i\p\j m 2) of the density 
matrix, 

Pmim.2 2k 

in which g rn denotes the "rate function" 



. _ ffroi + 9rn 2 

V 9rni + l9m 2 + lPmi+l,m 2 + l n Pmi,m 2 



(2.2) 



g m =j(j + l)-m(m-l). (2.3) 

The diagonal element p mm of the density matrix gives the probability to find the system of atoms in the state \jm); 
the elements p mi m 2 with mi — mi ^ will be referred to as coherences. It is worth noting a certain unidirectionality 
of the flow of probability and coherence, downwards the m-ladder, the physical o rigin of which is of course the 
low-temperature limit mentioned above. A further important feature of the system ( |2.2D is that the density matrix 
elements with different m\ — m-i evolve independently. To make that independence manifest it is convenient to 
introduce the quantum numbers 

mi+m 2 , mi - m 2 . . 

m = 2 ' 2 ( ' 

which can be simultaneously either integer or half-integer. Accounting for 

gmi + 2 9m2 =9m-k 2 (2.5) 

and changing the notation p mi m 2 f° r the density matrix element to pj^ we can rewrite the master equation as 

dp k 

-jf = 2k [y/g m+k +i9m-k+ip m+ i ~ {g m - k 2 )p k m ] . (2.6) 

It is now indeed obvious that the "skewness" k enters only as a parameter. 

The linear relation between the density matrices at the current time and at the initial moment, 

p k m (t) = J2 Dk mn(t)P k n(Q), (2-7) 
n 

defines the fc-dependent matrix D^ n (t) which will be called the dissipativ e pr opagator. Its column corresponding to 
a certain fixed n can be regarded as the solution of the master equation ( |2.6| ) corresponding to the initial condition 
Pm(0) — 5 mn . Due to the unidirectionality of the master equation it is obvious that flj,„ = if m > n. We shall 
drop the superscript k in the case k = 0, i.e. when the diagonal elements of the density matrix are considered. 
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III. SADDLE-POINT ASYMPTOTICS OF THE DISSIPATIVE PROPAGATOR 



There have been a number of successful attempts to treat the large- j limit of the superradiance problem ||,[) 10 7],||,[iJ ■ 
These were concerned with the solution of the master equation for certain particular cases or directly aimed at 
specific average properties of the process. The purpose of the present paper is to establish uniform asymptotics of 
the dissipative propagator without such restrictions. We use the exact solution of the master equation in the form of 
the Laplace integral which was obtained long ago Q but remained largely unexplored. Previously established results 
for the propagator, the distribution of delay times, and time dependent expectation values follow from our uniform 
asymptotic propagator. 

Before embarking on our proposed asymptotic adventure it is convenient to adopt the parameter 

v / 7CH r i)«i + 5 = J (3.1) 

as a measure of the "size" of the angular momentum; the semiclassical formulae to be established take a prettier form 
if we use J rather than j. 



A. Laplace representation of the exact propagator 

Following jP let us recall the Laplace integral representation of the propagator. Defining the Laplace image as 
T>m n (z) = J Q e~ zt D^ nn (t)dt we turn our master equation into a recursion relation with the easily found solution 



T> k (-A 1 TT y/9i-k9i+k . . 

2Ky/9m-k9m+k 2^ + 9l ~ « 

To get the dissipative propagator itself we invert the Laplace transform. Introducing a scaled time 

t = 2nJt (3.3) 

and the quantity 



(j + n)\(j — m)\ 

Wmn = j I yi = 

l=m+l 

we bring our propagator to the form 



n TT U +ny.{j - my. 

Qmn = \ \ 91 = J—, \TT- \7 (3-4) 

± A (7 + m (7 — n) 

l=m+l w > w 1 



d ^ {t) = y/Q m -k, n -kQ m+k , n+ k r b+to ° dv eTV/J -q — i : (3 5) 



/ \ V Q rrl ~k, n ~kQrn-\-k,n-\-k 

I, .. *" ' 

l=m 

where b should be larger than the largest pole in the denominator. 



B. Relation between densities and coherences 



An unexpected new result of the representation (3.5) is an identity connecting the propagators for the diagonal and 
for the off-diagonal elements of the density matrix, 



D k mn (r) = D mn (r) VQ™-k,n-kQ m+ k,n + k j ^ g) 



For the proof it is sufficient to shift the integration variable in (3.5) to v — v — k 2 . Alternatively, the connection 
between the diagonal and off-diagonal density matrix elements can be checked by entering the master equation with 
the ansatz 

k _ (j + m)! y/U-m-k)\(j-m + ky. j Kk 2 t _ 
Pm (J ra)\ V(j+m-fc)!( J+ m + fc)! e Pm[t) ' ( '- 7) 
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the new unknowns p m {t) then turn out to evolve in time like probabilities, i.e. to obey (2.6) for k — 0. 

The positive sign of the exponents in these relations between probabilities and coherences is not a misprint: the 
coherence p^, = p m +k,m-k does decay more slowly than the den sity = p m ,m- Moreover, there is no conflict with 
the nowadays popular phenomenon of accelerated decoherence |12| , |l3| , |l5| : Quantum dissipative processes do imply 
much larger decay rates for coherences than for probabilities but only so with respect to certain states which are 
distinguished by the process itself; for the dissipative process studied here such distinguished states are, for instance, 
coherent angular-momentum states [ pT|]To| but not the states \jm). 

A simple illustration of the statement just made may be helpful, even if it amounts to sidestepping to another 
dissipative process for an angular momentum, the one described by the master equation Q/9 = k{ [J z ,pJ z ] + [J z p, J z ]}. 
In that case the eigenstates \ jm) of J z are the distinguished ones as is obvious from p\ n = — 4k/c 2 pJ^ : The probabilities 
p® m are all conserved while the coherences have decay rates growing quadratically with the skewness k. 



C. Saddle-point evaluation of the Laplace integral 



The relation (3.6) between probabilities and coherences clearly allows us to confine the remaining investigation to 



the case k = 0, i.e. to the propagator of the densities. Our goal is to do the integral in the exact formula (3.5) in the 
limit of large J. To begin with, let us rewrite that formula for k = as 



D mn (r) = / e ^ n (v,r) dv (3 g) 



with the exponent 



Z mn (v,T) = tv/J- ^2 ln(v + gi) . (3.9) 

i— m 

Suppose now that n - m > 1. Then since the number of terms in the sum Z is proportional to n — m its value 
is generally also large, which fact suggests a saddle-point approximation. The stationary points of the exponent are 
given by the solutions for v of 

- 1 

Zmn=r/J-J2-—=Q. (3.10) 

l—m 

All roots of this saddle-point equation are real as is immediately seen by putting v = x + iy and separating the 



imaginary part. We further note that to the right of the largest pole u max = ™f< {— gi} of the integrand in (3.5) 



the sum in (3.1C) decreases monotonically from +oo to as v grows from w max to infinity. Therefore we have one 
and only one root vq in that domain. Its position depends on the time r: When r goes to zero t>o tends to infinity; 
conversely, for r — * oo the saddle point vq approaches the pole at v max . 
The second derivative with respect to v of the exponent, 

" 1 

Z mn = E [ v + 12 > ( 3 - H ) 

is positive for real v which means that the direction of steepest descent from the saddle is parallel to the imaginary 



axis. The saddle-point approximation for the integral (3.8) thus gives 

Qmn ^Z m „(v Q ,T) 



V ™mn 



D. Euler-Maclaurin estimates for the sums 



To render the expression (3.12) useful, we must evaluate the three sums in Z mn , Z mn , Z mn . The familiar Euler- 
Maclaurin summation formula f(k) « J f(x)dx + (f(m) + /(n))/2 comes to mind first but is not immediately 
suitable for our purpose. We rather employ a modified version which involves nothing but an integral; to compensate 
for the absence of the extra boundary terms the integration interval is extended, 
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n+l/2 



m— 1/2 



f(x)dx. 



(3.13) 



The accuracy o f bot h summation rules is the same for smooth summands f{k). 

In applying ( 3.13 ) to the sum in the saddle-point equation we rewrite the rate function as gi = J 2 — (I — 1/2) 2 , 
introduce the rescaled variables 



m — 1 



M = 



J 



and obtain 



n 1 

E — 

f-t v + gi 



n 



d.r 



Vvq + J 2 
J 



2Ja 



In 



i-l «o + ^ 2 — £ 2 
The saddle-point condition ( |3.10| ) thus takes the form 

1 (a + v)(a-fj,) 
t = — — In — - 

2a (a — v) (a + fi) 



(a + v)(a — fx) 



(a — v)(a + fi) 



(3.14) 



(3.15) 



(3.16) 



It determines a as a function of v, fi, and r. As already explained above, the single root of interest is positive and 
larger than the larger of \u\. 

Similary proceeding with the sums in Z mn and Z mn we find 



J 6 Z„ 



1 

2^ 



a 2 — /j 2 



a— a(/i,i/,r) 



Z mn (v, t) = J [r(a 2 - 1) - 2{v - fx) In J + 2{v - fi) - a(a, fi, v)] 
with the auxiliary function 

a(a, fi, v) = (y + a) ln(j/ + a) — (fi + a) ln(/i + a) 
— (a — v) ln(a — v) + (a — fi) ln(a — fi) . 



(3.17) 
(3.18) 

(3.19) 



We should comment on the slight asymmetry in the definitions of the macroscopic variables p. and v in ( 3.14 ). 
The use of (m — 1)/J instead of ra/J as the macro scopic variable /i is formally related to our extension by 1 of the 
integration interval in the summation formula (3.15) and has the benefit of preventing the small parameter 1/ J from 
appearing explicitly in the saddle-point equation (3.1 6|) . 



IV. UNIFORM ASYMPTOTICS OF THE PROPAGATOR 



We came to our saddle-point approximation assuming that the number of terms in the sum Z mn equal to n — m 
is large. It is not surprising therefore that the approximation (3.12) loses its accuracy when n — m is of the order 
unity or zero; that situation prevails, e.g., for small times r; an alternative approximation is then desirable and will 
be constructed presently. 



A. Small-time approximation 



To explain the essence of the new approximation let us give a simple example. Consider the Laplace image function 
with two simple poles V(z) = (z — c — d)~ 1 (z — c + d)" 1 and its original function V(t) = e ct dr 1 smhtd. As long as 
td <C 1 the hyperbolic sine can be replaced by its argument such that V(t) » te ct . We have thus in effect replaced 
the two close by poles of the Laplace image by a single second-order pole; that replacement is obviously justified for 
sufficiently small times. 

To employ this observation for the Laplace representation of the propagator ([T^) we introduce the new integration 
variable x — tv/J and obtain 



D,, 



i(t) = Q nm ( 



1 

2ri 



6+200 



"dx 



n?= m l* + 9iT/J] 



(4.1) 



G 



The length of the interval on which the poles of the integrand now lie is proportional to r, 

, r | m + n — 1 1 
J ~~ J 



(n 



(4.2) 



If that length is much smaller than unity the poles of the integrand of (4.1) are nearly degenerate, and that proximity 
enables us to replace the product in the denominator by the (n — m)-th power of the average factor x + gr / J with 

The integral is then easily calculated and yields the small-time asymptotics of the 



= J — 



re+m — 1 ^ 



9 = gap.-- V 2 
dissipative propagator, 



■Dmn ("7") 



Qr 



G)' 



exp ■ 



(4.3) 



Unlike the saddle-point approximation, the foregoing expression is fully explicit. We sha ll keep referring to it as the 
small-time approximant although the underlying small parameter is the combination (4.2) of both r and the quantum 
numbers m, n. 



B. Matching the two approximations 



The saddle-point and the small-time approximations for the propagator practically coincide for an intermediate 
range of arguments. Let us assume I = n — m+l^s> 1 but on the other hand £ = l/J -C 1 ( say, I ~ \/~J). The 
solution a(r,fj,, v) of the saddle-point equation (3.16) can then be found by expanding in powers of £, 



a 2 = v 2 + ^-^ + 0(C 2 ) 



(4.4) 



The exponent ( 3. IE ) in the saddle-point formula then simplifies according to r(a 2 -1)«( + t[(i/ - C/2) 2 — 1] and 
a(a, /z, v) sa £ (in | + 2 ) while the prefactor becomes S w r 2 /(. Collecting these pieces in (3.12) we obtain 



D„ 



V2^ ^l' W 



(4.5) 



l—n—m+l 



This in turn is the small-time approximation (4.3) provided we there replace the factorial (n — m)\ = (I — 1)! a la 
Stirling, (/ — 1)! w (5) • Hence the saddle-point and small-time approximations agree for 1 < ( C J. 



C. Uniform approximation 



The two app roxim ations under discussion can be merged into a single one which generally behaves like the saddle- 
point formula (3.12) but preserve s its accuracy even when m is close to n and/or the time r is small. We just have 
to divide the saddle-point result ( 3.12| ) by the ratio of the factorial (n — m)\ to its Stirling approximant. If n — m is 
l arge that ratio is unity but otherwise the correction replaces the saddle-point version with the small-time propagator 
(4.3). We thus obtain the principal result of our paper for the density propagator in the large-j limit, 



D„ 



QmnJ 3 ^ 2 ( le \ j[r(a 2 -l)- CT (q. M .i/)] 

(i-i)\VTe\J*J 

l = n — m + l, fj,= (m— 1)/J, v = n/J, a = a{fi,v,T). 



(4.6) 



It is valid in a wide range of quantum numbers and propagation times and thus merits the name uniformly asymptotic 
propagator. The error is of order l/J 2 except for the not very interesting late times when the bulk of the probability 
has settled in the lowest level; that latter restriction for t arises due to the close encounter of saddle and pole mentioned 
in Sect. IIIC. 



We have checked that (4.6) provides an efficient tool to numerically calculate the dissipative propagator; if j is large 
its accuracy becomes comparable or even superior to that of the numerical integration of the master equation. The 
only inconvenience is the necessity to determine the saddle-point parameter a = a{yL,v,r) by solving (3.16) which 
generally has to be done numerically. 
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V. SPECIAL CASES 



We proceed to considering situations in which the uniform approximation simplifies. The strategy invariably is 
to approximate factorials of large numbers a la Stirling. Some cases even allow for an analytical solution for the 
saddle-point parameter a whereupon fully explicit formulas for the propagator arise. Some well-known results of 
superradiance theory are thus recovered and revealed as special cases of the uniform approximation. 



A. Semiclassical approximation 



The uniformly asymptotic propagator (4.6) depends on the quantum numbers m,n 7 j in two ways. First there is 
the factorial dependence which reflects the discrete character of the representation. Second, there is the dependence 
on the arguments [A, v which can be regarded as the classical counterparts of m, n scaled with respect to the total 
angular momentum; they tend to continuous variables in the classical limit. 

Suppose we are not interested in effects tied up with the discreteness of quantum levels and want to obtain a smooth 
function of the macroscopic coordinates /i, v only. This is easily achieved by replacing the factorials (n — m)!, (j ± 
n)\, (j'±m)! by their Stirling estimates. While such a replacement would be unacceptably inaccurate if the arguments 
m,n approached ±j ("the poles" of the Bloch sphere in classical parlance) or each other, it otherwise reliably yields 

D m Jr) = 1 ; e J * ( "-"- r) , (5.1) 

$(//, v. t) = r(a 2 — 1) — a (a, fi, v) + cr(l, fi, v) . (5.2) 

We here speak of the semiclassical approximation because of the implied assumption that all the quantum numbers 
and their relevant combinations are large. As a function of [i at fixed v and r the semiclassical propagator displays a 
single maximum located according to 

6\X 1 — \JT 



i.e. a = 1. The saddle-point equation ( 3.16j ) then yields the most probable value of fi = J z /J at time r related to the 
initial value v through 

T=-ln-, 'r-) — r . (5.4) 

Written in terms of the polar angle of the Bloch vector co s 8 = fi, cos 00 = ^ the last equation becomes the solution 



of the equation of motion of the overdamped pendulum (1.1) mentioned in the Introduction. Indeed, the classical 



picture of the atomic dynamics in superradiance is that of the Bloch vector creeping from whatever initial orientation 
#o towards the equilibrium 8 — n like an overdamped pendulum with the azimutal angle <f> = arctan( J x / J y ) fixed. 
None too surprisingly, the maximum of the distribution D mn with respect to m occurs at the point m = m(n, r) 
predicted by the c lass ical motion of the Bloch vector. 

As it stands in ( |5.l| ) the semiclassical propagator correctly describes a broadening of the initially sharp distribution 
D mn (T — 0) = 6 mn to one with a width cx \/~J . For many applications that width is negligible such that we may 
replace the propagator by 

lim JD mn (r) = 5{pL - /i(r, v)) (5.5) 

J — >oc 



where Ii(t,v) is the classical trajectory according to (5.4). For instance, expectation values like (J^_J^Jf_} can be 



calculated to leading order in J with the help of the foregoing sharp version of the semiclassical propagator through 
the integrals 

( J s + J k z J-) - J 2s+k J % - Mr, "))(1 - l* 2 )y k dfx = (1 - m(t, v) 2 ) 3 ^ vf , (5.6) 

provided, it is well to repeat, the initial point n — Jv is well removed from the most highly excited ones, j — n> 1. 
No quantum effects at all survive in that expression; they would only show up as small standard deviations at most 



of order 1/y/J if the small width of the propagator (5.1) were kept 
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B. Early stage of superradiant decay of highest-energy initial states 



We now take up the previously best studied aspect of superradiance, the decay of the most highly excited atomic 
initial states, j — n -C j . We begin by studying the early stage, i.e. small t, while the bulk of the probability still 
resides with highly excited states. This means that only those propagator elements are significantly different from 
zero for which the final quantum number m is also close to j, or j — m <C j. 

We are so led to examine our uniform approximat ion w hen the macroscopic variables v and fj, are close to unity. 
Expanding the solution of the saddle-point equation ( 3.1 6| ) in powers of 1 — i>, 1 — /i we find the function a(/i, v, r) in 



terms of the nonlinearly rescaled time 

£ = e~ 2T (5.7) 



as a ps ■ From here it is easy to establish the ingredients of the uniform propagator (4.6), 
(a 2 - l)r - a ps (1 - v) ln£ + (y - /z) [ln(l - f) - ln(f - /z) - ln2e] , 

Hps- , Q m „ ps (2J)™ m ^ , (5.8) 

u-fi (j-ny. 

which bring the propagator to the limiting form 

D mn {r) = Ql^)e" n+1 (l - n ' m , (5-9) 
known as the linear approximation describing the early stages of the superradiant process O . 



C. Bright stage of superradiant decay of highly excited initial states 



Suppose now that the initial level is close to but the final quantum number m far away from j such that j — m 
is of the order of j. For simplicity we shall also assume that m is not close to — j. In classical terms, we take the 
Bloch vector as initially pointing almost to the north pole, but we wait long enough for it to develop a substantial 
component transverse to the polar orientation, i.e. a strong dipole moment; by excluding the late stages of near 
south polar orientation we confine ourselves to the phase of brightest radiation which actually gave rise to the term 
"super" radiance. 



Under the limitations on m, n just specified the saddle-point equation (3.16) can still be solved analytically. The 



important fact is that the function a takes on values close to unity. More accurately, it can be shown that the difference 
1 — a is of the same order of magnitude as 

S u = l-v, (5.10) 
the deviation of the initial classical coordinate from unity. It will be convenient to introduce the quantum time shift 

t' = t - T c i ass (n, v) , (5-11) 
where T c i ass (fi, v) denotes the classical time of travel from v to \i given by ( |5.4| ); in the situation under study it is 

/'- 



Tela 



2 5 V 2 1 



/' 



We can now write a as 



-It' 



S u . 



(5.12) 



(5.13) 



By similarly evaluating the other ingredients in the propagator (4.6) to leading order in 5 U and in addition replacing 
all factorials but (j — n)\ =11 by their Stirling estimates we come to 



D„ 



i+i 



J(l-M 2 ) 



I ! 



exp 



-2(Z + 1)t' 



-It' 



(5.14) 



1=3 -n 
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To connect with wellknown results we ban the quantum time shift r' by substituting J5.ll ), (5.12) and introduce the 
rescaled variables 



z = 2Jc 

The propagator thus assumes the equivalent form 

2 

Dmn (t) — 



x l+1 e x 



J(l-At 2 ) l\ 



(5.15) 



(5.16) 



l—j—n, fi=(rn—X)/J 



The special case of full initial excitation, I = j — n = 0, yields a distribution first derived by De Giorgio and Ghielmetti 

0* ... „ 

Contact with several previous treatments of superradiance is made by considering the bright-stage propagator ( 5.16| ) 
for high initial excitation as a function D(p, r; n) of the final coordinate \x and the time r and verifying it to obey the 
first-order partial differential equation 



dD 



d_ 

Oft 



(5.17) 



Obviously, that dynamics is devoid of quantum effects: The propagator D drifts along the characteristics of 5.17] , i.e. 
the fully classical trajectories 5.4 



D{h,t; n) 



1 -a* 2 



D(i/(fj,,T),0;n) 



(5.18) 



where ^(/x,r) is the time reversed classical trajectory obtained by solving (5.4) for v. All quantum effects inherent in 
the superradiant pulses then originate solely from an effective initial distribution D(fi, 0; n) which we read from (5.1£) 
by there setting t = 0, 



D((i,0;n) = 



J(l- M 2 )(j-n)! V 1 + M 



2 J 



l-/i 



j'-n+l 



exp — 2J 



1-M 
1 + M 



(5.19) 



We should emphasize that this effective initial distribution does not coincide with the true sharp initial form of 
the propagator, simply because our asymptotic propagator (5.16) is not valid at small times. The essence of the 
earlier theories of Refs. is thus recovered: Each run of a superradiant decay of a highly excited atomic initial 

state produces a macroscopic, i.e. classical radiation pulse originating from effectively random initial data, the latter 
reflecting quantum fluctuations. 



D. Time dependent expectation values 



We shall here establish a master formula for the set of "moments" defined as 

M fc .(r;0=tr [p(j - I; r) J» J*J S _] 



(5.20) 



with nonnegative integers k,s,l and p(j — 1;t) the density operator originating from the pure initial state — I). 
In the case of j much greater than 1 and k, s, I much smaller than j the average Mk s {T, I) can be written in the form 
of an integral over the classical variable \i with the propagator J£) mj j_;(r) = D(fi, u,t) as a weight, 



M ks {r-l) = J 



2s+k 



D( M)i ,,r)(l- M 2 )V^ 



(5.21) 



Upon employing the propagator (5.16) pertinent to the most highly excited initial states, changing the integration 
variable to x [cf. (5.15)], and once more using the rescaled time z from ( |5.15 ) we recover 



M ks (T-,l) 



J 



2s+k 



{4z) s e z 



x x l+s (z-xf 
o (z + x) 2s+k 



c dx . 



(5.22) 



an asymptotic result found by rather different methods in (|JT^]. It has a scaling form inasmuch as Mks(T;l)J 2s+k 
depends on J and r only through the single combination z. 
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VI. PASSAGE TIME DISTRIBUTION 



In the classical picture of superradiance the Bloch vector starts its downward motion from a certain initial angle 
©o and crosses the latitude at a strictly definite time r c / ass (/U, v) with v = cos©o, // = cos0. In other words, the 
classical probability density of the times of crossing a given coordinate p on the way fro m the initial point v is given 
by the delta function S(r — T c i ass (n, v)) = <5(r'); the quantum time shift r' defined in ( 5.11 ) is strictly zero in the 
classical limit. 

Let us now introduce the quant um mechanical generalization of the classically sharp passage time distribution. 
According to the master equation ( [2.6] ) for the densities, the change of the probability for the system to be in level m 
during the time interval dr equals {g m +iPm+i — 9mPm) dr. The quantity g m Pm{j)dT is obviously the probability for 
the atoms to go down from level m to level m — 1 during the time interval [r, r + dr] and 



Pm{r) = g m Pm(T) 



(6.1) 



is the corresponding probability density for the time of passage through level m. In particular, by stipulating the 
atoms to have started from the pure state \jn) with n > m, we specify the passage time distribution as proportional 
to the propagator, 

P m (r; n) = g m D„ m (r). (6.2) 

By simply integrating D mn (T) as given by the master equation ( |2.6| ) one easily shows that our passage time distribution 
is properly normalized to unity, 



P m (T;n)dT = 1. 



(6.3) 



Our uniform approximation for the propagator allows to easily and accurately calculate the passage time distri- 
bution. In particular, if the initial state is not close to the north pole, the function P m (r;n) is just a somewhat 
widened variant of the classical delta distribution, with a width inversely proportional to the square root of the 



second derivative J^^j. at the maximum of the exponent in the semiclassical approximation (5.1) 



However, for the more interesting initial states of highest excitation, the passage time distribution has little in 
common with its classical analogue. As follows from (5.14) in the case of the initial state \j,j — I) with I « j we 
rather get 



P m {T;j-l) = \(l + \ 



exp 



-2(1 + 1)t' 



'4 



-It' 



(6.4) 



This density depends only on I and r'. It gives directly the time distribution of the m — > m — 1 transition with respect 
to the classical time which corresponds to r' = 0. 

The absence of any explicit dependence on m and j means that the time distributions of probability calculated 
for different values of these quantum numbers but the same I = j — n differ only by a trivial time shift equal to 
the change in the classical time r c ; ass . In particular, the standard deviation of the time of crossing the mth level, 



At = V< 



< t > 2 with 



< T k >= 



T k P m {r-j - l)dr, 



(6.5) 



is a function of I only. The integrals (6.5) are easily calculated and give the mean passage time and the standard 
deviation as 



< t > 



Tela 



In I 



1 1 
^ k 

k=l 



At = 




(6.6) 



where C = 0.5772156649 ... is Euler's constant; in the case / = th e sums over k are absent. 

When I becomes large compared with unity the distribution (6.4) becomes sharply peaked around the point r' = 
predicted by the classical theory. However, as long as I remains of order unity or even becomes zero as for complete 
initial excitation the passage time distribution is rather broad: The relative standard deviation Ar/(r) is of order 
1/lnj; the small initial quantum uncertainty of the polarization sin (9 « a ^/Vj is found to be amplified to 
macroscopic magnitude in the passage time. 
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VII. APPENDIX: UNIFORM AND SEMICLASSICAL APPROXIMATIONS FOR THE PROPAGATOR 

OF COHERENCES 



The uniform approximation for the dissipative propagator of the non-diagonal elements (fc 7^ 0) is obtained via 
the exact relation ([0]). In the semiclassical approximation Stirling's formula is also applied in order to replace 
Qm±k, Qn±k by smooth functions of macroscopic arguments. We first note the uniform approximation for the prop- 
agator of coherences 

m — k , n —kQ rn + k , n + k / elY J 3 ' 

rnn ~ (J -1)1 Uv Vie 

xexp{j[r( a 2 -l + fc 2 /J 2 )- ( T( a , M , l /)]}| ;=n _ m+1 . (7.1) 

The saddle-point parameter a and the functions a, S do not depend on fc and are determined in exactly the same way 
as for the density propagator. 

Finally, we note the semiclassical approximation extending (5.1) to the propagation of coherences. Since there is 
an additional quantum number fc whose range goes to infinity when j — ► 00, a new macroscopic variable rj — k/ J has 
to be introduced. It is notationally convenient to write the previously incurred function cr(a, \i, v) with the help of a 
new auxiliary function 

q(x, y) = (x + y) \n{x +y)-(x-y) ln{x - y) ; (7.2) 

as cr{a,n, v) = q(a, v) — q(a 1 fi). Thus equipped we can present the propagator of the elements of the density matrix 
with skewness k as 



D k = 1 ■■ 



./<!>' 



run 



x /[l-/i-r / ) 2 ][l-( / i + 7 ? ) 2 ]V2^JS 



= - [q(l, v + 77) - q{l, fi + 77) + q{l, v - ry) - q(l, [x - r])] 

-a(a,n, v) +r(a 2 - 1 + if) . (7.3) 
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